Steady States of a Microwave Irradiated Quantum Hall Gas 
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We consider effects of a long-wavelength disorder potential on the Zero Conductance State (ZCS) 
of the microwave-irradiated 2D electron gas. Assuming a uniform Hall conductivity, we construct a 
Lyapunov functional and derive stability conditions on the domain structure of the photo-generated 
fields. We solve the resulting equations for a general one-dimensional and certain two-dimensional 
disorder potentials, and find non-zero conductances, photo- voltages, and circulating dissipative cur- 
rents. In contrast, weak white noise disorder does not destroy the ZCS, but induces mesoscopic 
current fluctuations. 
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The observation of giant magnetoresistance oscilla- 
tions in a microwave-irradiated two-dimensional electron 
gas (2DEG) 0], has spurred intensive theoretical activ- 
ity. Two distinct microscopic mechanisms for conduc- 
tivity corrections have been proposed: (i) The displace- 
ment photocurrent (DP)0, which is caused by photo- 
excitation of electrons into displaced guiding centers and 
(ii) the distribution function (DF) mechanism, which in- 
volves redistribution of intra-Landau level population for 
large inelastic lifetimes j3(|. 

Andreev, Aleiner and MillisQ have noted that ir- 
respective of microscopic details, once the radiation is 
strong enough to render the local conductivity negative, 
the system as a whole will break into domains of photo- 
generated fields and spontaneous Hall currents. In their 
proposed domain phase, motion of domain walls can ac- 
commodate the external voltage, resulting in a Zero Con- 
ductance State (ZCS) in the Corbino geometry, or a Zero 
Resistance State for the Hall bar geometry, in apparent 
agreement with experimental reports Ij. However, one 
may ask, what should be the effects of long-wavelength 
(relative to the cyclotron radius) disorder, which is ei- 
ther naturally present or deliberately introduced? What 
is the nature of the coupling between a disorder potential 
and the photogenerated fields , and could the disorder 
pin domain walls? Such pinning would affect the macro- 
scopic transport and could destroy the ZCS. 

In this Letter we incorporate a long-wavelength dis- 
order potential 0d(x) into the non- linear magneto- 
transport equations. We explore its effects on the domain 
structure and macroscopic transport coefficients. For the 
case of a constant Hall conductivity, we construct a Lya- 
punov functional (Oj which greatly simplifies the determi- 
nation of the stable steady states and their conductance. 
We use it to derive general stability conditions on do- 
main walls in the strong radiation regime. We also show 
that weak 'white- noise' disorder is an irrelevant perturba- 
tion, which does not destroy the ZCS. It docs introduce, 
however, mesoscopic non linear current fluctuations. We 
find solutions for the following disorder potentials: (i) 
The general one- dimensional potential, where domain 



walls are pinned to the potential extrema, which results 
in a non zero conductance and photo-voltage, (ii) The 
simple "egg-carton" potential solved variationally, and 
(iii) a generic non-separable potential depicted in Fig. ^ 
which is solved numerically, (ii) and (iii) exhibit two- 
dimensional domain wall pinning and frustration effects, 
which result in circulating dissipative currents. 




FIG. 1: Domain phase in the presence of disorder: the pho- 
togenerated potential 4>{x, y) corresponding to a generic long- 
wavelength disorder potential cj>d. The steady state is deter- 
mined by numerically minimizing the Lyapunov functional 
(see text). The domain walls (potential edges) are pinned by 
ipd, yielding a finite conductance. 



Non-Linear Magneto-Transport. In the presence of an 
external microwave field, we use a local relation between 
the dc current j(r) and the local electrostatic field E(r): 



j=j'^(E,r) + (7^zxE 



(1) 



Here we assume at the outset that the Hall conductivity 
is a constant, independent of r and E, which leads to 
considerable simplifications. The dissipative conductiv- 
(E) = dji/dEp, satisfies 



(2) 
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The vector function j'', in general, will depend explic- 
itly on the position r, due, e.g., to inhomogeneities in 
the 2DEG, and its direction may not be perfectly aligned 
with E. Eq. is supplemented by the continuity equa- 
tion, V • j = — p, where p is the charge density. We 
emphasize that equation contains all microscopic in- 
teractions at lengthscales shorter than the cyclotron ra- 
dius Ic, which serves as an ultra-violet cut-off, of order 
1/im. 

Writing E = —W(f>, we may relate changes in the elec- 
trostatic potential (p to changes in p through the inverse 
capacitance matrix W: 



S(t){r) = / d^r'W{r,r')6p{r'). 



(3) 



If a time-independent steady state is reached, then we 
have simply V • j = 0, and the precise form of W is 
unimportant. 

In a Corbino geometry, one specifies the potential on 
the inner and outer boundaries of the sample, and one 
looks for a solution for (j){r) consistent with these bound- 
ary conditions. Since we asssume to be a constant, 
the Hall current cannot contribute to V • j in the interior 
of the sample, so it does not appear in Kirchoff's equa- 
tions. Consequently, the solution for (j){r) is independent 
of and we may, for simplicity set = 0. To re- 
cover the Hall current, one inserts the solution E into 
the second term in (^. 

Condition ||2Jl on j'' allows us to define a scalar Lya- 
punov functional as 

G[<P] = J dVg(E), 9 = 1^^ dE'-f{E'). (4) 
A variation of l@J is given by 



5G = I (fr\/ -f 5(t)- [ dsn-f 

J J bou7id 



(5) 



The second integral vanishes on equipotential bound- 
aries, or in the absence of external currents. The extrcma 
of G are found to be steady states, with V • j = 0. Using 
the positivity of the inverse capacitance matrix W, one 
may show that G[4i{t)] is indeed a Lyapunov functional, 
i.e. a non-increasing function of time, so that its minima 
are stable steady states. In general, G may have multi- 
ple minima. Any initial choice of (f>{r) will relax to some 
local minimum of G, but not necessarily the "ground 
state" with lowest G. Nevertheless, we expect that in 
the presence of noise, the system might tend to escape 
from high-lying minima and wind up in a state with G 
close to the absolute minimum. 

Using the boundary term in |(SJ, the current across a 
Corbino sample is equal to the first derivative of G with 
respect to the potential difference V between two edges, 
and the differential conductance is given by the stiffness, 
or the second derivative: 



The Domain Phase. We now consider a homogeneous 
system, in the regime of strong microwave radiation at 
frequencies slightly larger than the cyclotron harmonics 
Lu > mujc,m — 1,2, .. . i.e. positive detuning. Both DP 
and DF mechanisms produce a regime of negative con- 
ductivity j(E) • E < 0, which implies a minimum of g{E) 
at a finite field |E| — Ec > 0, which was estimatcd j] 
to be of order hwc/{elc)- In order to satisfy equipoten- 
tial boundary conditions, and the constraints § dl-E — 0, 
field discontinuities and charge density singularities must 
form. 

A second order expansion of g about E^ reads as 

go(E) = go{E,) + ]^{E - E,fa, + A|V • Ej^ . (7) 

The clean system of Eq.©, is governed by a 'Mexican 
hat' Lyapunov density, with a fiat valley along |E| — Ec, 
i.e. the steady state local conductivity is 'marginally' sta- 
ble everywhere except inside the domain walls. The field- 
derivative coefficient A w '^Jdw implements the ultra- 
violet cut off, introducing a domain-wall thickness scale 
Idw assumed here to be of the order of 1^. Domain walls 
yield a positive contribution to G of order acE^ldw per 
unit length. In the absence of disorder, the system will 
simply minimize total domain walls contribution, subject 
to aforementioned constraints. 

A change in the average field (E) , required if there is a 
change in the applied voltage V, can be accommodated 
by a motion of domain walls, or a reorientation of the 
local E. The relative corrections to zero conductance 
vanish as Idw over the sample length. This defines the 
clean ZCS phase described in Ref. 

Long-Wavelength Disorder. In an inhomogeneous sys- 
tem, there will be a non-zero electrostatic field, Ed{r) = 
— V(/)(i(r), present in the thermal equilibrium state, with 
no microwave radiation or bias voltage. We may ask how 
this disorder field will alter the ZCS. 

At weak disorder field, jEdj << Ec, the Lyapunov den- 
sity near E Ec is modified to 



5(E, Ed) = 5o(E) - ai{E) E • Erf(r) 
which yields a current density 



0{Ej), 



nr) = -.,E, + t^:^E, 



(8) 



(9) 



dl 
dV 



fG 
dV^' 



(6) 



where X' = dX/dE, and the coefficient ci depends on 
microscopic mechanisms. 

We wish to elaborate on a physical issue regarding 
Eq.®: In the non-irradiated (dark) linear response the- 
ory, the current is driven by the electrochemical gradi- 
ent e = E — Erf. Similarly, one may expect the pho- 
tocurrent of the DF mechanism to also depend on e . In 
contrast, the 'upstream' photocurrent, pumped by the 
DP mechanism, involves transitions between single par- 
ticle states which feel the local electric field E(r). Thus, 
due to both contributions, even if the DP mechanism 
is relatively weak, E^ cannot be eliminated from Eq.Q 
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by a change of variables E ^ e . By our microscopic 
estimatejg, for the pure DP mechanism, CTc,cri(£'c) are 
close to the dark conductivity, and a'l « (Tc/Ec- 

Local stability requires that (T'^(r), of Eq. has non- 
negative eigenvalues. The lower (transverse) eigenvalue 
is given by 



■go - q-jE ■ Ec; 
E 



0{Edf > 



(10) 



so marginal stability occurs at E = Ec + <7[Eid ■ E / ac- 

In a steady state, the normal current density (in di- 
rection n) is continuous across a domain wall. If Ei and 
E2 are the fields on its two sides, we find by (|10|l and® 
that 



cr_(£;i) Ei-n = (7_(£;2) E2-n +0{Ej) 



(11) 



When Ei • n, i = 1,2 have opposite signs (as they do in 
the clean limit) Hll() can only be satisfied for a-{Ei) = 
a-{E2) = OiEdf' . This restricts the fields at the domain 
wall to be at their respective marginally stable values. As 
a result, the current density ij^ at the domain wall re- 
duces to 



j = -aiiE,)-Ed + OiEd 



(12) 



By Ea. (|12|l . current conservation and Gauss' theorem, 
and we obtain a global condition on any closed domain 
walls. 







nj 



(13) 



where Q^^ is the integral of the "2D disorder-charge den- 
sity," V • E(i/27r, over the area enclosed by the loop. 

Finally, we note that generically, the differential con- 
ductance of a sample in the Corbino geometry can be 
obtained by solving for the conductance of a linear sys- 
tem with local conductivity given by (T'*(E(r)), in series 
with resistive elements along the domain walls, which 
arise from movement of the domain walls in response to 
a variation in the applied bias V. (There could also be 
discontinuities in the current at discrete values of V , if 
the system jumps discontinuously from one local mini- 
mum of G to another.) We shall see that for weak long- 
wavelength disorder, the scale of the macroscopic con- 
ductance is set by the domain-wall contribution. 

White-Noise Disorder. In the ZCS, we now show that 
a weak 'white-noise' disorder potential, with a correlation 
length £^d (of the order Idw h), and root- mean-square 
value (t>d ^ Ec^d is an irrelevant perturbation which does 
not introduce new domain walls or destroy the ZCS. This 
is shown by using an Imry-Ma comparison^ of surface 
to bulk contributions to the Lyapunov functional. By 
for a square domain of area L^, the negative contri- 
bution of aligning E with the averaged disorder field E^, 
scales as —(Ji(j>dEcV Ed£,d- However the linear cost of its 
domain walls grows as +acEcLdldw Therefore weak dis- 
order cannot necessarily break the system up into smaller 
domains. 



The disorder, however, will produce current fluctua- 
tions across pre-existing domain walls, needed to satisfy 
boundary conditions on the sample. By Eq. (|12|) . the 
current density integrates across the domain wall to yield 
a random number of order 



SI = ±<Tl 



SVcpd ^1/2 
p ^5/2 



(14) 



While the conductance will average out to zero at volt- 
ages V > Ec^d or if multiple domains are in series. The 
random currents and conductance fluctuations should be 
observable in small samples, or as harmonic noise gener- 
ation for oscillatory bias voltage. 




-y J y 



FIG. 2: a) Photogenerated potential solutions for one di- 
mensional disorder <j)<i{y), with zero current and zero voltage 
boundary conditions, b) Map of V • E (red-positive, blue- 
negative) for solution (j>{x,y) of Fig. 1, showing the domain 
walls. Circulating dissipative currents are illustrated in one 
domain. 



One- Dimensional Disorder. In contrast to weak white- 
noise disorder, potentials with long-range correlations 
can be relevant perturbations. Consider the case of a gen- 
eral one-dimensional disorder (see FigEt), where (jidiy) is 
independent of the x-coordinate. At wavelengths larger 
than Idw , the Lyapunov functional is minimized if the 
system breaks up into parallel domains, so that E is ev- 
erywhere aligned with E^, and j'^ = 0. These conditions 
determine E(y) via (O, and the boundary conditions for 
a rectangular Corbino geometry (periodic boundary con- 
ditions at a; = and x = L). The voltage difference V 
between the leads at y = and y — W , satisfies 



V{I) 



w 



dy{E{y) - Ed{y)) 



(15) 



At strong radiation intensity and zero current, 11211 and 
transverse stability imposes that domain walls form pre- 
cisely at maxima and minima of (pd, given by the yi,i = 
1, . . . , TV To lowest order in (pd, there will be a non-zero 
photovoltage which only depends on these positions: 



^(0) = E, 




N 




0{Ed), (16) 
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where i = 1,3. . . are maxima. The differential conduc- 
tivity ayy is given, to lowest order in (jj^i, by 

1 _ 2Ec y-^ I _ 2Ec , , 

If P{Ed) is the probability distribution for Ed at a ran- 
dom point, it can be shown that E^^ = P(0). If 
Ed is taken from a Gaussian distribution, then Ed — 
(27r)^/^i5™'', where E™^ is the root-mean-square value 
of Ed- If 4>d is a single sine wave, then Ed = 2^/^7r_E™. 
The transverse differential conductivity Uxx can also be 
calculated, using H1U|I . and is given, to first order in (jjd 
by cr^x = cTi (|Erf|)/i?c. For a Gaussian distribution, 
one has <Jxx — (2/7r)crj,j,, while for a single sine wave, 

Uxx = (4/7r2)(Tyj^. 

Two- Dimensional Potentials The simplest 2D choice 
for (pd is the separable 'egg-carton' potential, given by 

E^ 

"Pd--^ (cos(x) + cos(2/)) . (18) 

We construct a zeroth order trial solution for zero bias 
current Jq = by placing domain walls on the lines x ~ 
WK and y = mn, for integer n and m. This yields constant 
electric fields in each square domain, of the form: 

Eo = ^(±£±y). (19) 

Eo is a gradient of a continuous potential, and satisfies 
charge neutrality H13|l . Upon application of an external 
voltage in an arbitrary direction, domain walls will move, 
as in the one-dimensional case, and also tilt with E into 
a herringbone pattern. A variational calculation, assum- 
ing that within each domain E is constant, finds, to first 
order in Ed, that axx = Q.83aiE^/ Ec- Numerical cal- 
culations confirm that the variational solution is at least 
close to the exact answer. 

We have calculated analytically the first order (in Ed) 
corrections to (|19|) . for zero external voltage, by integrat- 
ing Kirchoff's lawsf^. Away from the domain walls, we 



find E > Ec, which corresponds to circulating dissipa- 
tive currents, which match onto the tangential currents 
at domain walls, given by Eq. ^1 In Fig. ^ a generic 
two dimensional example is displayed, (jjd contains 20 
Fourier components chosen from Gaussian distributions 
with < |0rf(kp > independent of k. The potential 4> is 
found by numerically minimizing G. Fig. plots the 
2D charge density where domain walls appear as line sin- 
gularities. In both two dimensional examples, (|18|l and 
FigH G is frustrated from perfect alignment of E and 
Ed(r) by the conditions |E| > Ec and V x E = 0. This 
frustration underlies the circulating dissipative currents 
which are illustrated in one of the domains in Fig. 

In summary, we have introduced long wavelength dis- 
order into the transport theory of the microwave irradi- 
ated quantum Hall gas, using the Lyapunov functional as 
an organizing principle for the stability of steady states. 
We showed that weak white noise disorder is irrelevant 
for the stability of the ZCS although it produces meso- 
scopic current fiuctuations. For a strong and long range 
potentials, the ZCS state breaks up into domains, which 
will generally result in a finite conductance and a pho- 
tovoltaic effect. A microscopic theory necessary for the 
steady states dependence on microwave power and detun- 
ing frequency is deferred to a forthcoming publication^. 
We have not considered effects of conductivity anisotropy 
and variations in Hall coefficient. The latter will not af- 
fect the domain pattern or the longitudinal conductivity 
for one-dimensional disorder, but might have large effects 
and be experimentally relevant in other geometries. 
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